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Abstract: In this work we present a large area of aspects related to the problem of sampling 
inspection in the case of reliability. First we discuss the actual status of this domain, mentioning 
the newest approaches (from a technical view point) such as HALT and HASS and the statistical 
perspective. After a brief description of the general procedure in sampling inspection, we offer 
what we did call here as „personalized procedures": this means that we take into 
consideration the specific statistical law for time-to-failure. 

The original part refers to the construction of (n, 0) sampling plans and the use of ISO 
standard 2859 (MILSTD 105 E, the American original) in order to derive sampling plans by 
linking the AQL indicator (a defective fraction, in fact) to the hazard rate function. 

Some examples are given and necessary tables are provided also. 


Key words: reliability sampling plans; AQL; HALT; HASS; personalized procedures; (n,0) - 
plans; power distribution; Rayleigh law 


1. Preliminaries: the actual status of the problem 
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Reliability and its main component - durability - are considered as dynamical 
quality characteristics in the sense that possible nonconformities of the underlying entities 
(products, components, systems) are put into light in the process in which these ones perform 
their functions. For instance, we may observe the exploitation of the object by a client or may 
organize specifical lab durability/reliability tests. 
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A truck company for instance, records the mileage of its „moving units" up to their 
first failure and these records may constitute a database for the reliability analysis of that 
peculiar truck brand. 

Laboratory tests are more specific since the experimenter must choose a certain 
strategy for them: to use complete or (multi) censored samples, a time-truncation procedure 
(all performed under normal or accelerated conditions) a.s.o. 

The cost element must be involved too in such test designs. 

A huge monograph of Blischke and Murthy (2000, [3]^) describes new research 
directions in this field and provides a detailed list of ISO, NATO and from other professional 
associations (e. g. SAE - Society of Automotive Engineers, U.S.A.) documents regarding 
reliability and related domains (see also an earlier work of Kochendarfer and Pabst, Jr. 
(1971) [18]). 

It is important to mention that in the TQM (Total Quality Management) frame work 
or in the recent SIX SIGMA approach (see Gupta, 2004 [12]) reliability is often considered as 
a specific technical aspect since: 

(i) does not appear directly in the manufacturing process; 

(ii) is imbedded in the design stage of an entity in a provisional form; 

(iii) the acceptance (receiving inspection) of final products is usually performed 
on the so-called static characteristics such as length, volume, pressure, voltage, hardness 
a.s.o.; 

(iv) to test and to evaluate reliability requires an adequate logistics and a 
documented statistical know-how (that is a qualified personnel). 

It is also necessary to remember that most of the complex products generated by 
nowadays industries are „mecatronic entities" - that is their structure is a combination 
between mechanical, electrical and electronic parts each one with its own degree of 
automation and self-adjustment: this structural and functional complexity needs high 
engineering techniques for testing as well as an appropriate methodology for data analysis 
since the components of a mecatronic system have their specific behaviour in exploitation. 

In this sense, the classical techniques - such as running in of various vehicles at the 
highest level of stress intensity allowed by specifications - are combined now with two new 
procedures known as HALT and HASS. Proposed by dr. Gregg Hobbs in 1995 at QUALAAARK 
Corp. in Santa Clara (California) and registered as a federal mark these represent 
advanced methodologies for accelerated life-testing in order to determine as quick as 
possible the item reliability: this „quick evaluation" is needed to evaluate the possibility of a 
rapid intervention in the design stage, the goal being functional improvement (see also 
Dowling, 1999, [7]). 

HALT or Highly Accelerated Life Test means an experiment in which the objects are 
submitted to a highly accelerated test in order to estimate the durability of that kind of object 
that is its failure-free working period. 

The method is used in the design stage in order to detect weak points of the 
intended item to by mass produced, such as those probable nonconformities be eliminated 
before launching: the corrective action consists here in re-designing the product. 

HASS - or Highly Accelerated Stress Screening means an experiment in which an 
item is submitted under accelerated stress in the manufacturing stage in order to detect and 
eliminate the so-called „hidden defects", thus, preventing them to be genetically transmitted 
by the process to future batches. 
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The elements (stresses) used in HALT and HASS procedures are rapid transitions of 
temperature and omni axial vibrations (with six degrees of freedom) - this combination of 
random vibrations and thermo-cycles being the most performing way to accelerate the 
failure mechanisms of any kind. 

The idea to combine various types of stresses in not new: the engineers working in 
the „resistence des materiaux" (T. Albert, 1 829 and J. V. Poncelet, 1 839 in France, P. Hodge, 
1850 in England or A. Z. Wohlen and Z. Bauwesen, 1858, in Germany) have used this 
combination of random mechanical stresses with various temperatures in order to investigate 
the durability of metallic materials employed in railroad industry (details in Cioclov, 1975 [5] 
page 5-10). 

Some authors believe that this study of material strength face to various random 
mechanical and thermo stresses can be considered as the embryo of the reliability theory 
(see Barsan-Pipu et al, 1999, [2]). 

A last remark in this context: accelerated tests are needed especially if we deal with 
very high reliable objects for which is almost impossible to wait their natural failure in order 
to evaluate numerically the underlying reliability. The exclusive usage of complete samples 
and normal conditions has proved to be in many cases time consuming, extremely expensive 
or even not very significant from a technical viewpoint (examples of such pitfalls are given in 
Meeker and Hamada, 1995, [20]). 

As regarads the receiving inspection of products for which the characteristics of 
interest is durability/reliability, we distinguish several approaches which are quite different 
since are based on distinct principles. These are: 

a) Attributive approach: viewed as the oldest and simplest one, it considers the 
durability/reliability as a measurable characteristic which can be attributivisable (an 
attribute, as for instance conforming or nonconforming). In our framework, we submit to a 
reliability test a witness batch of size Ng (non-repairable items, for example), during a fixed 
period of time (Tg); then, record the number (d) of failed elements over that period and 
compute the estimated defective fraction p = d/Ng of the lot. This value (p) may by taken 
as the desired AQL needed by the ISO standard 2859 or its American variant MIL STD 105 E 
in order to employ for instance single sampling plans (n, A | R = A -I- 1), where A and R are 

acceptance and rejection numbers (see Iliescu, 1982, [13]). 

The major advantage of attributive method is its simplicity: one has to compare the 
number (d) of failed items in a sample of size (n) - this one given by the standard, with the 
acceptance number A. If d < A the lot is accepted - otherwise, that isd>R = A-l-l, the lot 
is rejected (this is the case of single sampling plans). 

The main disadvantage in this case is the reductionism itself of the method: it does 
not take into account the specific law for failure behaviour, which is a key element in 
reliability evaluation (Catuneanu - Mihalache, 1990, [4]). 

As a consequence, the so-called personalized procedures have been devised and 
these acceptance schemes use effectively the statistical distribution of time-to-failure. 

b) Average operating time approach: this procedure takes into consideration a 
specific time-to-failure law and establishes acceptable and unacceptable values for mean 
operating time. Using a sample of size (n) and an acceptable number A (which may be a 
minimum member of hours of failure-free operation, considered acceptable), the lot will be 
accepted or rejected in comparison with this number. 


192 




Quantitative Methods Inquires 


JOURNAL 

OF 

APPLIEO 

QUANTITATIVE 

METHODS 


c) Hazard rate approach: the procedure is similar to the above one, the fixed 
values being of hazard (failure) rate. 

d) Sequential approach: we have in mind Abraham WALD'S procedure which 
can be applied in both cases - attributive and variable ones (see Wald, 1973, [23]). 

In this paper we shall propose another procedure, namely linking the fraction 
defective (p) with the hazard rate, obviously taking into account the specific form of time-to- 
failure distribution. 


2. The general procedure 


Attribute sampling inspection may be performed regardless the very nature of 
quality characteristic tested: static or dynamical one. There are considered two proportions P, 
and P 2 of product units, first acceptable with 1 — a = 0,95 probability (that isa = 0,05) and 
the second one acceptable with a lot smaller probability, usually P = 0,10. We may have 
alsoa = p. The problem consists in determining the sampling plan (n. A) where n is the 

sample size and A is the acceptance number. To find n and A we proceed as follows: 

(i) from the relationship (see Baron et al, 1988, [1]). 

2 

^2 _ ^ \-P\m . . 

P ~ r" 

1 A a;m 

where is the £ - quantile of the chi-square distribution with m 

degrees of freedom (see [13]), we approximate „the best m" for which (1) is fulfilled. 


(ii) 

A = 


the acceptance number A is given by 
where [ ] is the integer part 




(iii) 

n = 


X 


the sample size (n) is furnished as 

2 

(rounded to the nearest integer) 


( 2 ) 

(3) 


Example: if we take P, = 0,25% and P 2 = 1,75% and a = P = 0,05 , based on 

chi-square tables, we get 


A 

1 

2 

3 

4 

5 

6 

7 

8 

P 2 /P 1 

13 

7,5 

5,7 

4,6 

4,0 

3,6 

3,3 

3,1 


Since P 2 /P 1 = 1 ,75/0,25 = 7, hence the nearest value is 7,5 for which we extract A = 
2. If m = 2 (A -I- 1) = 6, we shall obtain 

1.635 


n = 


Zo.05-6 


= 327 (units) 


(4) 


2x0.0025 0.005 

Therefore, the sampling plan is (n = 327, A = 2) which seems to be not very 
economical (since n is quite large) for an expensive testing. 


193 




JAQM 




JOURNAL 

OF 

APPLIEO 

QUANTITATIVE 

METHODS 


Quantitative Methods Inquires 


3. Personalized procedures 

We shall consider now the time-to-failure model, namely the mathematical object 

(5) 


'f(t;0), te[0,+oo), 0 = ( 0 j, 02 ,...eJ, eR, j = l,k 

00 

we havef(t; 0>O, Vt>0), jf(t;0)dt = l. 

0 

Here, T is the continuous random variable which represents the time-to-failure 
behaviour. The simplest model is considered; 

3.1. The exponential model 

This is described by the density 

T ; f (t; 0) = (l/0)exp(- 1/0), t > 0, 0 > 0 (6) 

where 0 has the significance of average durability, since E(t) = 0 and its inverse is just the 
hazard rate h(t;0) = f(t;0)/R(t;0) where R(t;0) = 1 — p(t;0) = exp(- 1/0) 
is the reliability function. 

In this case we shall fix two values for E(T), namely 0[ (acceptable mean-life with 1 — a 
associated probability, 0<a<l) and 0j (undesirable mean-life with P associated 
probability 0 < P < 1) where usually a = 0,05 and P = 0,10 . 

If we will employ an „r out of n" durability test, that is we will obtain a censored sample 

r<n (7) 

the average life estimate is 


1 

=- 

r 


Zt,+(n-r)t„ 


( 8 ) 


(see [1], vol. I, page 532). 

The acceptance numbers are as below 


Ai = 


01 -Xa;2r 


2r 


, if we adopt producer's risk variant or 


Aj — 


02 ■ Xl-p;2r 

2r 


, if we prefer customer's risk variant. 


The decision to accept the lot is taken if 0 ^ or 0^ — . 

Example; Let 0j = 1000 (hours), a = 0,05 and r — 4 when the sample size in 
n = 26 units (we work using producer's variant). From chi-square tables ([!]) we find 
Zoos-s ~ 2.733 and hence Aj = lOOOx 2.733/8 « 342 . 

Suppose that the average mean-life was 0^ =960 (hours). SinceO^ > A^, the 
underlying lot is accepted with confidence level of 95%. 
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3.2. The Weibull model 

Proposed in 1951 (see [24]) by a Swedish military engineer, Wallodi WEIBULL 
(1887 - 1979), this model is considered as a generalization of exponential and Rayleigh 
laws (see Isaic-Maniu, 1983, [14]). Its density is 

T;f(t;0,k) = (k/0)t^ ‘exp(-tVe} t > 0 , k > 0 (9) 

and has the following average value 
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E(t) = 0'* • r(l + l/k) where from 0 = 


E(T) 


r(i + i/k) 


l/k 


( 10 ) 


uu 

Here r(u) = | e * • t“ ’dt is the well-known Gamma function (see Dorin et al. 


1994 [6], page 239 -244). 

If p = Po is the proportion of nonconforming objects, then 


Po =l-exp^- 


t-r(l + l/k)' 

E(t) _ 


( 11 ) 


where 0 has been replaced by its expression from (10). 
Hence, we deduce 

t _ [-ln(l-po)]''^ 


( 12 ) 


E(t) E(1 + l/k) 

and taking t = Tq as testing time (expressed usually in hours) we can compute To/E(T) if k and 
Po are known. 

As it has been shown in [22, page 112] this approach leads often to large sample 
sizes, which is not always convenient. 

This handicap could be eliminated by constructing (n, 0) type sampling plans, that is 
plans where the acceptance number is zero. 

3.3. The Gamma model 

A random variable T has a Gamma density function if 


T:f(l;e.k) = ^ 




0 




exp(-t'‘/0), t > 0, 0,k>O 


(13) 


It is also a generalization of the exponential one (for k = 1), we get 

f(t;0,l) = (l/0)exp(-t/0), t > 0, 0>O). 

The average life is in this caseE(T) = k0 . If the shape parameter k is known, then 
if we fix an acceptable mean-value k0Q and a testing time Tq, such that with a given 
probability P, the average durability be al least k0g, using the ratio To/k0(, and P, one may 

deduce the needed sample size to perform the test (details in Gupta and Groll, 1961 [11] or 
Voda, 1 981 [22], page 116-119). 

Example; Suppose that we give A = 0 (the acceptance number) P = 0.95, k = 
2 and the average-life k0g =10.000 (hours) - that is0g =5.000. If Tq (testing time) is 

1.000 hours, thenTg/kOg =0.10. From the below table (reproduced from [22] page 119) 
we detect n = 1 70 units. Therefore 
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a) a sample of n = 170 item are submitted to a test over the period of 
To = 1 .000 hours; 

b) if there are no failures during this testing period, the lot is accepted; 

c) if there exists at least one failure, the lot is rejected, since the acceptance number 


IS zero. 


To/k0o 

4 

1.0 

0.05 

0.10 

0.05 

0 

4 

10 

170 

639 

1 

6 

16 

269 

1013 

2 

9 

22 

358 

1444 

3 

1 1 

27 

411 

1655 


4. (n, 0) - type sampling plans 

These plans, no matter which are the other input elements are invariantly based on 
the zero acceptance criterion, that is always A = 0. 

A way to construct such plans is to fix in advance the following elements; testing 
time Tq, lower bound for reliability R(t) evaluated for t = Tg or a lower bound for the average 
lifetime E(T) and the consumer risk P . 

This procedure has been described in a series of research reports of SVUSS (Stdtny 
V^kumny Ustav pro Stavbu Stroju/Bechovice, CSSR - State Research Institute for Machine 
Construction/ Bechovice, former Czechoslovakya - see Drimlovd (1970 [8], 1973 [10]) and 
Drimlovd and Zaludovd (1971) [9]) and it was applied for the exponential law. In this case a 

lower bound for the reliability R(t;0) = t>0, 0>O may be easily obtained from a 

limit fixed for the mean durability, since E(t) = 0 . 

The sample size is given as 

-Inp 


n = 


exp(-To/0o) 


(14) 


where 0g is the acceptable value for E(T). 
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Example: Take Tq = 100 (hours) and R|o„er (100) = 0.90 and P = 0.10 (or 

10%). Since - lnP« 2.302, we get easily n » 2.3/0.10 = 23 units. Hence, the sampling 
plan is (n = 23, A = 0). 

For the Weibull law, taking account the results from § 3.2. and § 3.3., we get 

-Inp 


n = 


-|k 


• r(i+i/k) 


(15) 


where E(T)io„er is the mean lifetime which we wish to accept with P probability. If the sample 
exhibits a smaller value for E(T) then E(T)io„er/ the lot is rejected with (l — P) probability (in 
this model, the shape parameter k is assumed to be known). 

Example: Let us consider a Rayleigh distribution (which is a Weibull one for k = 

2 - see for details the present authors 1998 [16]). We wish to reject batches with a mean 





Quantitative Methods Inquires 


JOURNAL 

OF 

APPLIEO 

QUANTITATIVE 

METHODS 


durability less than 2000 (hours) with a 90% probability (that is p = 0.10). The testing time 
was fixed as Tq = 500 (hours) and obviously, A = 0. 

In the next table, we present some values of (n) for various P and To/E(T), using 
formula (1 5). 


p 

Tn / E ( T ) 

0.1 

0.2 

0.4 

0.8 

0.05 

382 

95 

24 

6 

0.10 

293 

73 

18 

5 


In our case, Tq/E(T) = 500/2000 = 0.4 and we have for P = 0.10 the sample size n 
= 18. units. Therefore the sampling plan is (n = 18, A = 0). 

For the power distribution, namely 


/ 1 V 

T;F(t;5,b)= - , 0 < t < b, 5>0 
lb J 


(16) 


which is a generalization of the uniform one (which is obtained for 5 = 1) or a peculiar case 
of Sedrakian's one (see Sedrakian, 1968 [21] or the present authors, 1995 [15]); 


T;F(t;b,c,5,k) = l- 


1 - 




b-c 


(17) 


where 0<c<t<b, 5,k>0 (the power form is recovered if c = 0 and k = 1), we have 
immediately for (1 6): 


E(t) = — and Var(T) = - ^ 

^ ^ 5 + 1 ^ ^ (5 + l)^(5 + 2) 

Hence, if 5 is known, if we fix a lower acceptable bound for E(T), we get 

1 + 5 


' lower 


■E(T) 


lower 


and taking into account the results from [22 page 1 14], we finally obtain 

Inp 


n = 


In 


1 - 


5A 
5 + 1 


(18) 


(19) 


( 20 ) 


where A = To/E(t)i^„^^ 

In the below table we offer some values for n given some A and 

p = 0.01;0.05;0.10and5 = 3/2;2. 


a 

P 

A = T „/ E ( T ),_, 

0.2 

0.4 

0.6 

0.8 

3/2 

0.01 

54 

18 

9 

5 

0.05 

70 

23 

12 

7 

0.10 

108 

36 

18 

11 

2 

0.01 

128 

31 

13 

6 

0.05 

167 

40 

17 

8 

0.10 

256 

62 

26 

13 
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Example; Assume that 5 = 2 and we with to reject batches with a mean 
lifetime less than 3000 (hours) with a 99% probability (hence P = 0.01). The testing time Tg 
is 600 (hours). 

We have A = 600/3000 = 0.2 and for P = 0.01 and 5 = 2 we read in the table 
n = 128. Therefore the plan is (n = 128, A = 0). 

5. The use of MILSTD 105 E 

As we said in § 1 , the attributive sampling in practice makes use of the well-known 
document MILSTD 105 E or ISO variant ISOSTD 2859 (see Iliescu, 1982, [13, page 162 - 
180]). 

This standard does not refer to purely reliability elements such as testing time, 
failure/hazard rate, MTBF (Mean Time Between Failures), EOT (Effective Operating Time) 
a.s.o. (for specific English acronyms used in reliability theory, see Kovalenko, 1975 [19, page 
437 - 466]). The items are simply divided into two classes: conforming and 

nonconforming/defective ones - no matter what indicator is considered. 

In this paragraph we shall construct sampling plans by linking the lot defective 
fraction (p = AQL) with the specific hazard rate of those objects, using some input elements 
of MILSTD 105 E such as lot size (N), code-letter (CL) which will lead to the sample size (n). 


Our approach assumes that the failure model is of the type 


f (t; 0) = ^3'(t; 0)exp[- (p[t, 0)] 

(21) 

where ^3(t;0) > 0 for every t > 0 and 0 > 0 and 


00 

jf(t;0)dt = l 

0 

(22) 

The form (21) provides the distribution function 

F(t;0) = l-exp[-^3(t;0)] 

(23) 

the defective fraction p being just p = 1 — e . 


Taking logarithms, we have 

-ln(l-p) = ^(t;0) 

(24) 


and since the hazard rate is h(t;0) = f(t;0)/R(t;0) 

where R(1;0) = 1 — F(1;0) = ex p[^(t;0)], we obtain h(t;0) = ^'(t;0) and therefore (24) 
becomes 

-^3'(t;0)ln(l-p) = h(t;0)^3 (t;0) (25) 

If we take now - ^(t;0) = 0t^, then h (t; 0) = 2 0 1, t > 0 0 > 0 which is just the 

classical Rayleigh hazard rate (see Barsan-Pipu et al, 1999 [2, page 91)). 

In this case (25) has the form 

-21n(l-p) = t-h(t;0) 

Consider now the simplest situation when we know N (the lot size), Tg (the testing 
time - assumed to be the life span of the items) and the acceptable hazard rate h(t;0) for t 
= Tg, expressed in failures per hour. 
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In order to ease the computations, we present some values of lOOt h(t;0) linked 
with some preferential AQL values listed in MILSTD 105 E (table II A) - see also 
Kirkpatrick, 1970 [17, page 363]). 


Table (.) 


AQL 

(%) 

0.15 

0.25 

0.40 

0.65 

1.0 

1.5 

2.5 

4.0 

lOOt-h(t;0) 

0.30 

0.50 

0.80 

1.30 

2.01 

3.02 

5.6 

8.16 


In these instances, the sampling procedure is the following; 

1) Knowing N (lot size - let's say N = 930 units) and using the general inspection 
level II (suggested in most of the cases by the standard - see Table I „Sample size and code 
letters", [17, page 362]), we find the code letter CL = J; we shall denote IL as inspection 
level. 

2) from Table II A „Single sampling plans for normal inspection" (already 
mentioned) we have to draw from the lot a sample of size n = 80 units which will be 
submitted to the test. 

3) assuming that the testing time is To = 1000 (hours) and the acceptable hazard 
rate forTo is h (1000) = 0.000003 failures/hour, we deduce 

100 TO h (TO) = 100 xlOOO x 0.000003 - 3.0 

and from the table (.) we see that the nearest value is 3.02 wich indicates AQL = 

1 .5%. 

4) with CL = J and AQL = 1 .5%, from the same table II A we read the acceptance 
number A = 3 (consequently, the rejection one is R = A + 1 =4). 

Therefore, during the testing period Tq = 1000 (hours) we shall observe the number 
(d) of failed elements from the sample (n = 80 units). If d< A = 3, the lot is accepted - 
otherwise (that is d > R = 4) the lot is rejected. See the below scheme: 
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